home *** CD-ROM | disk | FTP | other *** search
/ SGI Hot Mix 17 / Hot Mix 17.iso / HM17_SGI / research / lib / obsolete / chi_sqr.pro < prev    next >
Text File  |  1997-07-08  |  1KB  |  61 lines
  1. ; $Id: chi_sqr.pro,v 1.2 1997/01/15 04:02:19 ali Exp $
  2. ;
  3. ;  Copyright (c) 1991-1997, Research Systems Inc.  All rights
  4. ;  reserved. Unauthorized reproduction prohibited.
  5.  
  6.  
  7. function chi_sqr, a,DF 
  8. ;+
  9. ; NAME:
  10. ;    CHI_SQR
  11. ;
  12. ; PURPOSE: 
  13. ;    CHI_SQR returns the cutoff value v such that
  14. ;
  15. ;        Probability(X > v) = a
  16. ;
  17. ;    where X is a random variable from the chi_sqr distribution with 
  18. ;    v degrees of freedom. 
  19. ;
  20. ; CATEGORY:
  21. ;    Statistics.
  22. ; 
  23. ; CALLING SEQUENCE:
  24. ;    Result = CHI_SQR(A, DF)
  25. ;
  26. ; INPUT:
  27. ;    A:    probability
  28. ;    DF:    degrees of freedom
  29. ;
  30. ; OUTPUT:
  31. ;    If a is between 0 and 1, then the cutoff value is returned.
  32. ;    Otherwise, -1 is returned to indicate an error.
  33. ;-
  34. if (a gt 1 or a lt 0) then return,-1  
  35. if a eq 0 THEN return, 1.e12
  36. if a eq 1 THEn return,0
  37.  
  38. if DF lt 0 THEN BEGIN
  39.   print,'chi_sqr-- degrees of freedom must be nonegative'
  40.   return, -1
  41. ENDIF
  42. case 1 of
  43.  DF EQ 1: UP=300.0 
  44.  DF EQ 2: UP =100.0 
  45.  DF GT 2 and  DF LE 5: UP=30.0
  46.  DF GT 5 and DF LE 14: UP = 20.0
  47.  ELSE: UP = 12.0
  48. ENDCASE
  49.  
  50. Below=0
  51.  
  52.  while chi_sqr1(UP,DF) LT 1- a DO BEGIN
  53.       Below = UP
  54.       UP = 2*UP
  55.  ENDWHILE
  56.  
  57.  
  58. return, pd_bisection([1-a,DF],'chi_sqr1',Up,Below)
  59. end
  60.  
  61.